MODULE EDDVAR
!
! Verwendungszweck:
!
!      Bereitstellung des Eddingtonfaktors und verwandter Variablen
!
!-----------------------------------------------------------------------
!
!     f_edd    : Eddington Faktor
!     mu_r     : Abstrahlungsrandbedingung radial (Standardwert 1/2, gilt fuer f_edd = 1/3 )
!     mu_z     : Abstrahlungsrandbedingung in z-Richtung (Standardwert 1/2, gilt fuer f_edd = 1/3 )
!
!-----------------------------------------------------------------------

      use config, only : rkind, np

      implicit none

      private
      public :: EDDFCT, f_edd, mu_r, mu_z


      real(kind=rkind) :: f_edd(np)
      real(kind=rkind) :: mu_r, mu_z

CONTAINS





SUBROUTINE EDDFCT
!
! Verwendungszweck:
!     Berechnung der Eddington-Faktoren fuer jede Zelle
!
!     HIER:  verwenden Eddington-Approximation:  f_edd = 1/3
!
!       wenn man spaeter einmal den eddingtonfaktor wirklich ausrechnet muss man
!       sich auch dessen zeitliche Abstufung ueberlegen:
!            - zeitliche zentrierung: f_edd, f_eddZ
!            - nach jedem iterationsschritt eddington faktoren updaten - oder erst nach konvergenz ?
!              ersteres ist eleganter (und stabiler ?) und erlaubt groessere zeitschritte, zweiteres ist schneller
!              - so wie Materialfunktionen ?!
!

      use physco, only : z13, z12, z0
      use zvar,   only : pseudo2D

      implicit none


      f_edd(:) = z13

      mu_r   = z12

      if (pseudo2D) then
         mu_z = z12
      else
         mu_z = z0
      end if

END SUBROUTINE EDDFCT


END MODULE EDDVAR
